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ABSTRACT 


'  When  an  asperity  traverses  at  a  relatively  high  speed  over  the  surface 
..of  a  medium,  the  resulting  friction  force  in  the  contact  zone  may  cause  a 
"hot  spot"  that  would  lead  eventually  to  "heat-checking"  or  thermomechanical 
cracking  in  the  medium.  The  asperity  traversing  speed,  the  load  distri¬ 
bution  of  the  asperity  over  the  contact  area  and  the  relative  dimensions  of 
the  contact  area  all  affect  the  stress  state  in  the  medium,  thus  the  initi¬ 
ation  of  fracture  therein.  Moreover,  the  mechanical  and  the  thermal  pro¬ 
perties  of  the  medium  also  play  significant  roles  in  the  susceptibility  of 
themomechanical  cracking  )b#-the~medi urn,' >thus  wear  of  the  friction  surface. 

In  the  present  paper,  it  Is  shown  that^by  assuming  a  uniform  distribution  of 
the  load  over  the  contact  zone,  the  predicted  stress  level  could  be  40%  less 
than  that  of  a  non-uniform  distribution  of  an  equal  total  load.  While  the 
shape  of  the  contact  area  does  not  affect  the  stress  level,  its  size  does. 
For  the  same  average  pressure,  increasing  the  size  in  the  direction  of 
asperity  travel  raises  the  stress  level;  while  increasing  the  aspect  ratio 
(relatively  increasing  the  length  normal  to  the  traversing  direction)  lowers 
the  stress  level.  In  the  limit  of  an  infinite  aspect  ratio  (a  plane  strain 
model),  less  than  20%  of  the  stress  level  at  the  aspect  ratio  of  one 
results.  The  multifacet  influence  of  asperity  speed  manifests  in  rate  of 
thermal  energy  input  and  in  the  thermal  layer  thickness,  which  determines 
the  magnitude  of  thermal  stress  and  the  location  of  fracture  initiation.'  'In 
the  mechanical  and  the  thermal  properties  of  the  medium  material ,  higher 
coefficient  of  Coulomb  friction  and  higher  Young's  modulus,  as  expected, 
lead  to  higher  stress.  But  higher  Poisson's  ratio  does  not  signif icantly 
Increase  the  stress.  High  thermal  expansion  coefficient  and  low  thermal 
capacity  are  well  known  to  be  the  source  of  high  thermal  stress.  The  result 


bears  out  their  effect  in  thermomechanical  cracking.  However,  thermal 
capacity  and  thermal  conductivity  jointly  Influence  the  thermal  field 
through  the  thermal  diffuslvity,  which  occurs  In  the  governing  differential 
equation.  The  report  demonstrates  that  the  thermal  diffuslvity  Is  Indeed  a 
derived  thermal  parameter.  It  Is  the  values  of  thermal  capacity  and  the 
thermal  conductivity  individually  that  determines  the  thermal  field,  thus 
the  stress  field  that  leads  to  thermomechanical  cracking  of  the  wear  medium. 
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NOMENCLATURE 


Asperity  characteristic  dimension,  the  half  width  of  the  rec¬ 
tangular  contact  area  In  the  direction  of  traverse  or  radius  of 
the  circular  contact  area. 

Half  length  of  the  rectangular  contact  area  perpendicular  to  the 
direction  of  traverse. 

Specific  heat 

Young's  Modulus 

Thermal  conductivity 

)  Load  distribution 

Average  pressure  over  the  contact  area 

Peclet  Number  (=  Va/k) 

Heat  flux  through  the  contact  area 

Traction  over  the  contact  area  in  the  x,  direction 

Aspect  ratio  (b/a)  or  time  coordinate 

Temperature  field 

Displacement  field 

Traverse  speed  of  asperity  (-x^  direction) 

Coordinates  fixed  to  the  moving  asperity 
Partial  derivative  with  respect  to  x^  coordinate 
Coefficient  of  thermal  expansion 
Dimensionless  temperature  field  (=  Tk/qQa) 

Thermal  diffusivlty 
Lame  coefficient 

Lame  coefficient,  modulus  of  rigidity 
Coulomb  coefficient  of  friction 
Poisson's  ratio 
Mass  density 


Stress  field 


M 

Mechanical  Stress  field 
<j|j  Thermal  stress  field 

aj  Maximum  combined  thermomechanical  principal  stress 

Oimensionless  coordinates  (=  x^/a) 

Dimensionless  depth  of  maximum  thermal  gradient 
(“)  Fourier  transform  of  (  ) 


1.0  INTRODUCTION 


In  some  mechanical  systems,  mating  surfaces  are  pressed  against  each 
other  and  undergo  relative  sliding  motion,  which  may  occur  with  dry  friction 
either  by  design  or  Inadvertently  through  normal  operation.  Had  the  pres¬ 
sure  been  evenly  distributed  according  to  design,  the  service  life  of  the 
device  would  not  be  a  serious  problem  even  at  a  high  rubbing  speed.  How¬ 
ever,  it  is  well  known  that  the  actual  contact  area  may  only  be  a  small 
fraction  of  the  nominal  area  at  the  design  interface  [1],  Burton  [2]  con¬ 
sidered  10'4  as  a  possible  areal  ratio  (the  ratio  between  the  actual  contact 
area  and  the  nominal  area  at  the  design  interface).  In  other  words,  a  low 
nominal  design  pressure  may  result  In  a  very  high  Interfacial  pressure,  thus 
a  very  high  dry  frictional  force  In  the  actual  contact  area.  Experiments 
have  demonstrated  such  high  friction,  and  the  resulting  extremely  high  local 
temperature,  called  a  "hot  spot"  [3,4].  The  hot  spot  is  considered  one  of 
the  causes  of  frictional  cracking,  or  "heat  checking"  [5].  It  suffices  to 
assume  that  the  localization  of  contact  area  is  due  to  some  form  of  asperi¬ 
ties,  from  thermoelastic  instability  or  impurities,  which  may  be  fixed  to 
any  of  the  mating  surfaces  or  may  precess  with  respect  to  both.  Based  on  a 
model  of  asperities  fixed  to  one  mating  surface  and  tneir  excitation, 
thermo-mechanical  actions,  on  the  other  hard  mating  surface,  Ju  and  Huang 
[6,7]  developed  a  three  dimensional  theory  of  thermo-mechanical  cracking. 

In  the  general  theory,  the  complicated  solutions  were  left  in  the  expres¬ 
sions  of  Fourier  transform.  Numerical  solutions  were  obtained  for  material 
properties  similar  to  those  of  Stellite  III.  It  was  shown  that,  for  a  con- 
servative  areal  ratio  of  10  ,  a  24d  kPa  (35  psi }  design  pressure  could 
result  in  a  240  MPa  (35000  psi)  asperity  pressure.  At  this  pressure,  a 
moderate  Coulomb  coefficient  of  0.5  and  an  asperity  traversing  speed  of  10 


mps  (400  ips)  would  lead  to  thermo-mechanical  cracking,  initiating  at  a 
depth  of  one-tenth  asperity  size,  approximately  100  ^m  deep.  A  general 
review  of  related  works  is  referred  to  [7].  The  present  report  shall  study 
the  asperity  and  the  material  parameters  which  Influence  the  thermo¬ 
mechanical  stress  state  in  the  medium  which  lead  to  heat  checking. 

1.1  Mathematical  Model  and  Parameters 

The  phenomenon  of  thermomechanical  cracking,  as  observed  from  experi¬ 
ments  and  operational  damages,  is  seen  to  associate  with  relatively  hard 
materials;  such  as  cast  iron,  Stellite  III  and  the  likes.  The  plastic  wear 
and  surface  shear  are  restricted  In  a  very  thin  surface  layer,  Blau  [9]  and 
Ruff  and  Blau  [10].  Thermo-mechanical  cracking  could  very  well  initiate  at 
sub-surface  region,  where  the  material  behavior  is  essentially  thermo¬ 
elastic.  The  governing  differential  equation  are  the  thermoelastic  Navier's 
equation  and  the  uncoupled  Fourier  equation,  respectively 

(\+{i)3^.Uj  +  =  (3\+2p)a9.jT,  (1) 

and 

k  T  -  T,  (2) 

where  3.  =  3/ax^,  (\,y)  are  the  Lame's  coefficients,  p  is  the  mass  density, 
a  is  the  coefficient  of  thermal  expansion  and  tr  the  thermal  diffusivlty.  The 
indices  in  the  subscript  i,j,k,  =  1,2,3  and  the  summation  convention  are 
used  for  all  repeated  indices  of  roman  minuscules.  Tne  stress  field  {a.^} 
can  be  expressed  in  terms  of  the  thermoelastic  Hookian  Law: 

=  \3|<U|<5ij  +  u(  3  j  +  3^  Uj)  -  (3\  +  2y )  aT  5  ^  j  ,  (3) 

where  4^  is  the  Kronecker  delta.  The  coordinates,  { x ^ ,  in  Equations 
(1-3),  are  material  coordinates;  that  is,  they  are  fixed  to  the  wear 


material  over  which  the  asperity  traverses.  The  field  variables,  {u^}  and 
T,  as  evidenced  in  Equations  (1,2),  are  defined  in  (x.,t).  In  other  words, 
time  is  an  explicity  variable.  The  moving  boundary  conditions  are  also  time 
dependent.  We  can  now  restate  our  mathematical  model  as  follows.  The  size 
of  the  asperities  under  consideration  are  of  the  order  of  1  mm;  the  thick¬ 
ness  of  the  hard  material  that  was  acted  on  is  at  least  an  order  of  magni¬ 
tude  larger.  Hence,  the  material  is  adequately  represented  by  a  half  space 
with  the  asperity  traversing  over  the  surface  boundary  at  a  uniform  speed 
(V)  as  shown  in  Figure  1.  At  high  speed  traversing,  the  high  temperature 
and  surface  yield  due  to  the  excitation  are  subgranular.  If  the  region  of 
consideration  is  away  from  such  subgranular  layer,  the  half  space  region  is 
basically  elastic  and  subjected  to  brittle  fracture.  For  the  present  pur¬ 
pose,  the  material  can  be  considered  as  homogeneous  and  Isotropic  without 
local  flaws.  In  the  analytical  formulation,  with  respect  to  the  material 
reference  frame  (fixed  to  the  half  space  material),  the  field  variables 
(u^,T)  will  have  local  time  variation.  Hence,  the  problem  is  transient. 
However,  by  taking  into  consideration  the  uniform  properties  of  the 
materials,  asperity  motions  and  loadings,  we  can  justify  the  Invariant 
states  of  the  field  variables.  In  other  words,  with  respect  to  a  convective 
reference  frame  (fixed  to  the  asperity  thus  moving  relative  to  the 
material),  the  analytical  formulation  becomes  "steady  state";  that  is,  there 
is  no  explicit  time  variable.  The  coordinates  shown  in  Figure  1  are  con¬ 
vective  such  that  X|  is  oriented  along  and  opposite  to  the  asperity- 
traversing  direction;  x^  is  perpendicular  to  the  boundary  surface  into  the 
material;  and  is  parallel  to  the  surface  but  perpendicular  to  the  tra¬ 
versing  direction.  The  governing  differential  equations  for  the  displace¬ 
ment  f^ld  (u^  and  the  temperature  field  (T)  in  the  material  referring  to 
the  convective  coordinates,  the  dynamic  thermoelastic  Navler's  equation  (1) 
and  the  uncoupled  Fourier  law  (2),  become: 
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(1)  OCTANGULAR  CONTACT 
AREA 


2a 


(\  +  +  ]i3jjU^  =  pV  +  ( 3\  +  2}i)a3jT,  (4) 

k3ijT  =  VijT.  (5) 

The  Equations  (4,5)  are  only  implicitly  time  dependent.  The  boundary  con¬ 
ditions  in  terms  of  the  convective  coordinates  are  also  time-independent.  On 
the  surface 

a3i  =  "  Ri’ 
and 

k33T  =  -  VRp  (7) 

where  the  tractions  are:  =  0,  R3  =  P(XpX2)  and  =  ^Pfx^Xo)  in  the 

contact  region  and  zero  everywhere,  P^pX^)  is  the  asperity  pressure, 

tif  is  the  coefficient  of  Coulomb  friction  and  k  is  the  thermal  conductivity. 

Without  loss  of  generality,  the  mating  surface  and  the  non-contact  wear 

★ 

surface  are  assumed  to  be  adiabatic.  The  regularity  conditions  for  the 
stress  and  temperature  fields  hold  at  Infinity;  that  is,  at  x^  =  ®, 

°-fj  =  (8) 

and 

T  =  0.  (9) 

From  the  governing  differential  equations  (4,5)  and  the  boundary  conditions. 
Equations  (6,7),  we  can  thus  delineate  the  parameters  as  those  character¬ 
izing  the  asperity  and  those  the  material.  The  groups  of  parameters  are 
listed  as  follows 


*  The  thermal  energy  loss  across  the  free  surface  by  convection  is  small 
generally.  Since  the  maximum  temperature  gradient  occurs  at  the  immediate 
vicinity  of  the  moving  asperity,  the  heat  loss  by  convection  will  nbt 
alter  significantly  on  the  criterion  for  thermo-mechanical  cracking. 
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Table  I 


V 

P(xrx2) 

a 

b 


Asperity  parameters 
Asperity  traversing  speed 
Asperity  pressure 

Asperity  contact  area  width  (In  the  traversing 
direction) 

Contact  area  length 
Contact  area  shape 


\,  } x  or  £,v 
p 

a 

if 

k 


Table  II 

Material  parameters 
Mechanical  constitutive  coefficients 
Mass  density 

Coefficient  of  Coulomb  friction 
Coefficient  of  thermal  expansion 
Thermal  dlffuslvlty 
Thermal  conductivity 


The  present  report  will  study  the  effect  of  these  parameters. 
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1.2  Analytical  Solutions 

For  general  solutions  of  Equations  (A, 5),  Fourier  transforms  of  co¬ 
ordinates  ( x ^ x 2 )  are  used,  resulting  in  ordinary  differential  equations  in 
x3;  that  is 

00 

7(x1,x2,x3)  =5f[f(x1)]  =  ^Jf  f(x1)exp[1{x1x1+  x2x2)]dx1dx2,  (8) 


where,  and  in  the  sequel,  (  )denote$  the  Fourier  transform  ?[(  )].  The 
temperature  fields  are  obtained  in  the  transformed  space  and  are  employed 
directly  as  input  to  the  transformed  differential  equation  of  Equation  (4). 
The  Inverse  transform  of  special  cases  are  carried  out  by  numerical  methods 
The  numerical  method  itself  Is  checked  with  the  solutions  of  temperature 
field,  which  can  be  alternatively  obtained  by  the  Green's  function  solution 
[8]. 

The  non-homogeneous  differential  equation  (4),  subjected  to  nonhomo- 
geneous  boundary  condition  (6),  is  solved  by  dividing  the  stress  and  the 
displacement  fields  {cr^,  u^ }  Into  two  portions 

|j;i,  u,  }  and  {a^  ,  u!}  such  that 

ij  1  ij  ' 

M  T 

^  =  u1  +  ui  ,  (9) 

°1j  =  aijM  +  a1jT’  •  (10) 

The  individual  displacement  and  stress  components  are  related  following  the 
thermoelastic  Hookian  Law,  Equation  (3),  as: 


"l/  *  xVk\;j  +  i*  (Vi”  +  V/)-  (: 

o-i -  X3kuk^s  1  j  +  11  (*jui^  +  ( 

M  M 

The  set  ,  u^  }  is  required  to  satisfy  the  homogeneous  differential 
equation  but  nonhomgeneous  boundary  condition;  that  is,  in  the  region, 

(\  +  +  u3jjU^  =  p V  3yjU.j  ,  (. 

and  on  the  boundary, 


°3i  =  ~Ri »  on  surface; 


a,  i  =0.  at  x,.x..  = 


The  set  u^}  Is  to  satisfy  the  non-homogeneous  differential  equation 

but  homogeneous  boundary  condition;  that  is,  In  the  region 

(\  +  n]9-jjUj^"  +  liSjjU^  =  pV^3|^U^  +  (3\  +  Zji)a3^T,  (15) 

and  on  the  boundary, 

a.^  =  0,  on  the  surface;! 

>  (16) 

a^T  =  0,  at  xkxk  =  ®.  ) 

The  mathematical  delineation,  for  convenience  of  solution,  essentially  in 
the  physical  sense  separates  the  stress  and  the  displacement  fields  into  two 
parts.  One  results  from  the  mechanical  excitation  of  the 

asperity;  those  are  the  traversing  normal  pressure  and  frictional  traction. 
Another  (<r^T,  u^}  is  solely  the  consequence  of  temperature  field  or  rather 
its  gradient  (3.T)X  resulted  from  the  heat  input  due  to  frictional  energy 
loss  rate.  They  are  respectively  called  the  mechanical  portion  and  the 
thermal  portion  of  the  solutions.  It  is  noticed  from  the  Fourier  equation 
(5,7),  that  the  temperature  field  is  itself  the  result  of  the  friction  force 
^P.  Friction,  therefore,  directly  affects  the  mechanical  portion  as  well 
as  the  thermal  portion  of  the  solutions.  The  relative  dominance  of  each 
portion,  as  evidence  from  the  thermal  boundary  condition  Equation  (7), 
depends  on  the  traversing  speed  V  of  the  asperity. 

The  general  solutions  of  the  mechanical  portion  in  the  Fourier  trans¬ 
formed  expressions  as  given  in  [6]  are: 


T»T'.; 


2jiF  u^M  =  (Ix^n^  -  -  x^r^l^)  e 


n  -n,x^ 

-C-  *  ~  -  K  \  ~  1  J 


{'*lnln2*3  ■ 


x-  XX  - 

-  [  £  <"3  -  2nl"2)  +  "2n3  ]  R1  +  ("3  -  2nin2>  ^  e 


—  M  — •  —  —  —9  ill  ao  /  ^ 

2]iF  =  (^2n3^3  *  *l*2n2^1  "  *|n2^  e  1  ~  |^2nln2^3  + 

XX 

+  -TT  ("3  -  2"l"2>  *1  -  [  n*  <n3  ’  2nln2>  +  n2n3  ]  M  e 


2^F  u3M  =  (n-^n^  +  Ix^n^E}  +  1^2nln2^2 )  e 


"nlx3 


-  [n}(^  +  x2^)  ^3  +  ^in3^i  +  ^x2n3^2^  e 


-"2*3 


.  A  *---- 


id 


[•] 


=  8  [(x^  +  x2  -  ri|  )  +  x1t]  e  { x-| tn1  n9R^  + 


1  12  3 


+  1X|  [x22N2  +  n2n2]  ^  -  ix^x2N2R2}  e  2  3, 

=  ff  [(x^  +  x22  *  n^2)  K  +  x22]  e  1  3  -  {x22r^n2E3  - 


-  1xnx22N2^  +  1x2  [X|2N2  +  n2n3]  ^  e 


~n2x3 


=  5  [(x^  +  x22  -  r^2)  K  -  n^2]  e  3  3  + 

—  7  —  ?  —  —  —  ~n^x-5 

+  n2  [(x^  +  x2  )  n^3  +  in3  (x^  +  x2R2)]  e  L  , 

-<V3  _  _  _  1  _ 

=  x^x2B  e  -  {X|X2n^n2R3  +  £  x2  fn2n3  “ 


-  (x32  -  x22)  N2]  j  * 1  [n2n3  +  "  *2 *)  N2^  e 


-2  -  2 

1/  ^  w  ^  ' 


=  -Ix-^B  e  3  3  +  {ix2n3  (x^2  +  x22  +  n22)  $3 


-  2X|X2n^n2R3  -  [x^2N^  +  n3  (n22+  x22)  R2y  e  2  3/2, 

.f|  x 

s  -ix3n^  e  3  3  +  {ix^  (x^2  +  x,,2  +  n22)  $3  - 
[x2  +  n3  (n2  +  X|  ) ]  R^  -  2x^x2n^n2^2}  0  2  /25 

=  ( Vs-  Mr2712)1/2J 

=  x  x  _  i  m  2-  2 

xsxs  2  n2  1  ’ 

=  V/cr,  =  (\  +  2j*)/p,  c22  =  v/pf 


-  nln2xsxs’ 


K  =  \/2]i  =  v/(l  -  2v) 


I  =  +  ^n2*s^s’  r,s  =  1*2, 

= 

The  expressions  of  the  mechanical  portion  are  quite  general,  applicable  for 
any  traction  {R^}.  Specifically  for  the  problem  of  asperity  excitation,  R3 
Is  the  asperity  pressure  P(XpX2),  Rj  is  the  asperity  friction  and 
R2  *  0  In  the  contact  area  of  the  surface  and  zero  everywhere  else  on  the 
surface.  By  taking  Rj  and  R3  separately,  the  computation  can  be  much 
simplified  and  also  can  afford  a  comparison  of  stress  fields  as  effected  by 
normal  pressure  and  frictional  force. 

The  temperature  field  in  the  Fourier  transformed  expression  is: 


-u)X^ 

T  -  e  {[C^  cos(ex3)  +  C2  sin(ex3)]  + 
+  [C3  cos(ex3)  +  C4  sin(ex3)]  F2h 


where 


c  =  ~2m2  +  in 

3  2u(w2  +  02) 


r  =  20  ±  in 

2  20(<o2  +  02)’ 

r  =  -n  -  12o>2 

3  2<u(<u2  +  02)’ 

r  -  -n  +  1202 
-  a  o  > 

4  20(<u  +  0^) 


(s4  *  - 


-[ 


[l  t  - 


(19) 


,  V  ■ > '  vv 
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2-2-2  - 
s  -  Xj  +  *2 »  n  =  -Vx^/r, 

and  =  VTT^  /  k . 

The  temperature  fields  for  individual  cases  are  left  in  the  Fourier  trans¬ 
formed  expression  as  the  excitation  to  the  transformed  differential  equation 
of  the  thermal  portion  {<r^,  u^)-  of  the  thermoelastic  Navier's  equation 
(15,16).  The  temperature  fields,  in  inverse  transformed  expressions,  are 
used  as  a  double  check  of  the  numerical  procedure  against  the  Green's  func¬ 
tion  solutions  and  also  serve  to  correlate  the  temperature  gradient  with  the 
maximum  tensile  stress. 

The  expressions  of  the  thermal  component  {<j^,  u^)  are  very  complex, 
that  it  is  rendered  necessary  to  resort  to  numerical  methods  in  inverse 
Fourier  transforms  for  individual  values  of  parameters.  The  Fourier  trans¬ 
formed  expressions  of  thermal  components  are  given  as 


where 


Gg  =  -1 (wG  -  0H) ,  Hg  s  -1 («H  -  0G), 

al  =  a4  +  2  s  2(1  -  v)/(l  -  2v) , 
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a2  =  a4  +  1  s  1/(1  -  2v), 

~  (3a4  +  2)a  =  2a(l  +  v)/(l  -  2v), 
a4  =  \/y  =  2v/(l  -  2v). 


2a3a2s  u>9 


?  9 

[a^se  H  +  (s  -  a^sw  +  a^0fc)  G]. 


and  h2  =  t"s9  H  -  (s2  -  s<,  +  a^2)  G]. 

The  Fourier  transformed  expressions  of  the  thermal  stress  field  are: 
*n  r  _wXq  -sxn  -| 

=  1  [  (-b|H  +  b2G)  e  3  +  (-b3H  +  b4G)  e  3  J  , 


-  T 
a22  . 


[  WAn  -i 

(-b5H  +  b6G)  e  3  +  (-b7H  +  bgG)  e  J  , 


-  T 
*33  =  i 


("bgH  +  b1(JG)  e  3  +  (bnH  -  b^G)  e  3  J  , 


-  T 

*12  _  - 


[UiA  *5  O  A  O  ** 

(-^23^  +  b^G)  e  +  ("b^H  +  bj^G)  e  J  , 


-  T 
q23 

-  >"  -  —  —  _  V/ 


[WA**  -■ 

(byH  -  b18G)  e  -  (b^gH  -  b2gG)  e  J  , 


-  T 

*31  _  - 


[  IDAn  -JAn  -■ 

(b^H  -  b|gG)  e  -  (bjgH  -  b2gG)  e  J  , 


*vi  v.  V  ^  V* 


lLklL  ^U»kLl>iL' 


where 


-  2  . 


>1  =  a^e  f2u,e  COS(0X3)  -  xi  sin(ex3)]s 


3  rr  2 


b2  =  a^?  txl  COS(0X3)  +  2“9  sin(ex3)j. 

b  =  !i  /!l  +  !iL .  xi2x3\ 

3  u  \&2^  aia2  al  '  ' 


a3  r  a|“  “  s  _  2  a4(s<u  “  s  ) 
4  a^(o8  |_  a2S  1  82 


+  (s  -  »)  x22x3  j  , 
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_3  rr  2 
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•m -- 
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M  M 

Equations  (17,18,20,21)  are  the  general  solutions  of  {a^  ,  u^  )  and 

u/}  in  their  Fourier  tranformed  expressions.  With  given  excitation 
conditions:  {R^,V}.  The  inverse  transform  can  be  carried  out  numerically. 
The  numerical  technique  was  reported  by  Ju  and  Huang  in  [6].  The  effects  of 
Individual  parameters  are  therefore  analyzed  through  specific  numerical 
variations  of  their  values. 


I 


2.0  ASPERITY  PARAMETERS 
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The  asperity  parameters  involve  those  excitation-related,  the  actual 
pressure  and  its  traversing  speed  (P,V),  and  those  related  to  the  geometry 
of  the  contact  area,  its  width  in  the  direction  of  traverse,  it  length  and 
its  shape  (a,  b,  contact  are  shape).  The  excitation-related  parameters  are 
principally  design  oriented.  Yet  the  actual  pressure  and  its  distribution 
in  the  contact  area  have  to  be  determined  empirically.  The  contact  area 
configuration  is  not  readily  known  to  be  controllable  by  design.  However, 
studies  of  these  parameters  are  important  in  analysis  of  wear  resistance,  to 
promote  reliabililty  in  design.  The  studies  will  provide  the  unstanding  in 
adequacy  of  analysis  by  using  a  two-dimensional  model  or  by  assuming  an 
uniform  distribution  of  pressure  {thus  friction)  in  the  contact  region. 

2.1  Excitation  Parameters 

The  traversing  speed  of  the  asperity  (V)  influences  both  the  convective 
acceleration  term,  Equations  (13,15),  and  the  thermal  input,  Equation  (7). 

It  is  readily  conceivable  that  at  low  rubbing  speed  the  mechanical  portion 
of  the  stress  dominates.  In  the  limit  at  zero  traversing  speed,  there  may 
still  be  mechanical  field;  but  the  thermal  field  becomes  trivial.  At  high 
speed  asperity  excitation  the  thermal  field  prevails.  The  present  report 
addresses  specifically  the  latter  case. 

At  high  speed  exultation,  the  temperature  is  high  near  the  surface,  and 
the  thermal  stress  is  compressive  as  a  result  of  constraint  against  thermal 
expansion.  At  below  surface  layer,  where  the  temperature  field  has  high 
gradient,  the  interparticle  constraint  thus  creates  high  tensile  stresses, 
which  at  unfavorable  loading  could  lead  to  crack  initiation  through  cohesive 
failure.  The  phenomenon  is  illustrated  in  Figure  2,  showing  the  stress 


field  near  the  surface  layer  as  a  result  of  high  speed  asperity  traversing 
at  15  m/s  (600  ips).  The  excessively  high  compressive  stress  at  the  surface 
layer  indicates  a  plastic  zone  up  to  a  depth  of  the  order  of  10  pm.  The 
depth  of  plastic  deformation  is  well  of  the  same  order  for  hard  wear 
materials  [9].  However,  if  the  crack  initiation  due  to  cohesive  failure 
occurs  at  a  depth  of  an  order  of  magnitude  larger,  the  thermoelastic  formu¬ 
lation  of  the  thermomechanical  cracking  model  in  Art.  1.1  will  hold.  The 
maximum  tensile  stress  due  to  thermal  field  will  occur  in  the  neighborhood 
of  the  depth  of  maximum  temperature  gradient.  The  temperature  field  and  its 
gradient  are  respectively: 


4>(£,n,S)  =  /  f (C,n, a )  exp(-c2/29)de,  (22) 

■'0 

=  -C  f *  9"1  f(5,n,9)  exp(-c2/29)d9,  (23) 

where  (£,n,t)  =  (x^/a)  is  the  dimensionless  coordinate,  <j>  =  Tk/qQa  Is  the 
dimensionless  temperature  field  In  the  medium,  qQ  is  the  heat  flux  due  to 
the  average  frictional  load,  and 


In  the  expression  of  f(£,n,9),  N  is  the  number  of  repeated  asperity  exci¬ 
tations,  t  is  the  aspect  ratio,  £r  *  Z  -  (r  -  1)\‘  -  1/2  Pe  9,  Is  the 
periodic  length  Vx/a,  x  is  the  period  of  excitation,  Pe  =  Va/r  is  the  Peclet 
number.  The  expression  is,  without  loss  of  generality,  using  a  uniform  load 


distribution  over  a  rectangular  contact  area.  The  depth  (c  ),  at  which  the 
maximum  temperature  gradient  occurs,  can  be  obtained  by  extremizing  the 
temperature  gradient,  Equation  (23),  with  respect  to  the  depth  c,  such  that 

9-2  (e2  -  9)  f(e,n,9)  exp(-c2/29 )de  =  0,  (24) 

S=<m 

The  maximum  temperature-gradient  occurs  at  a  depth  depending  on  the  Peclet 
number,  Pe,  which  is  controlled  by  the  asperity  traversing  speed  (V)  and  the 
thermal  diffusivity  {<).  The  solution  to  Equation  (24)  are  obtained  numer¬ 
ically  for  different  values  of  number  of  asperity  repetition  (N),  lateral 
position  (n) ,  contact  area  aspect  ratio  (t),  and  Peclet  number  (Pe). 

Figure  3  illustrates  the  depth  of  maximum  temperature  gradient  in  a  hard 
wear  material  such  as  Stellite  III  as  a  function  of  the  Peclet  number  for  a 
single  asperity,  n  =  0  and  an  aspect  ratio  of  10.  The  plot  is  up  to  some 
Peclet  number  corresponding  to  the  magnitude  of  the  traversing  speed  (V)  in 
the  neighborhood  of  the  Rayleigh  wave  speed  (cR),  beyond  which  the  mathe¬ 
matical  model  no  longer  holds.  It  is  noticed  that  up  to  Pe  =  1.4  x  10^  the 
depth  is  fairly  constant  around  0.1.  Based  on  the  Stellite  III  as  the 
medium,  an  asperity  width  of  0.5  mm  (-0.02  in.)  and  a  traversing  speed  of 
15  m/s  (-600  ips),  result  In  Pe  =  1200.  The  depth  thus  defines  the 
thermal  layer  near  the  wear  surface.  Since  the  layer  Is  indeed  by  order  of 
magnitude  larger  than  the  plastic  layer,  the  mathematical  model  is  justi¬ 
fied. 

The  effect  of  the  asperity  traversing  speed  may  thus  be  summarized  as 
follows: 

1.  higher  rate  of  frictional  energy  input,  resulting  in  higher  temper¬ 
ature,  thus  higher  thermal  stress, 


.100  CASE  I  (rectangular,  uniform  pressure) 


maximum  temperature  gradient  vs  Peclet  no. 


2.  reduced  thickness  in  thermal  layer  at  higher  speed,  leading  toward 
higher  gradient  thus  higher  thermal  stress. 

The  dominance  of  thermal  stress  over  the  mechanical  stress  at  high¬ 
speed  asperity  excitation  (15  m/s)  is  shown  in  Figure  4.  For  a  circular 
contact  area  with  paraboloidal  distribution  of  pressure  and  friction, 


the  thermal  principal  stress  is  six  times  higher  than  the  mechanical 
principal  stress.  With  uniform  distribution  of  pressure  the  ratio  is  almost 
five  times.  The  significance  of  the  thermal  components  of  the  stress  field 
cannot  be  overemphasized. 

The  load  distribution  in  the  contact  area,  taking  Into  consideration  of 
the  thermo-mechanical  behavior  of  the  wear  material,  has  not  been  determined 
either  analytically  or  experimentally.  With  the  result  of  Hertzian  contact 
problem,  the  assumption  of  uniform  distribution  of  loading  Is  indeed  an 
oversimplification.  The  analysis  assumed  a  paraboloidal  distribution  with 
the  same  total  load  as  the  uniformly  distributed  load.  Equation  (22).  The 
difference  in  the  thermal  principal  stresses  of  the  two  is  already  evident 
in  Figure  4.  Figure  5  further  compares  the  total  thermomechanical  principal 
stress  of  a  paraboloidal  distribution  and  a  uniform  distribution  of  loading 
in  a  circular  contact  area.  For  the  same  total  asperity  loading,  the 
assumption  of  a  paraboloidal  distribution  would  result  in  a  thermo¬ 
mechanical  principal  stress  around  80£  higher  than  that  resulted  from  the 
assumption  of  uniform  distribution.  Of  course  neither  may  claim  to  be  the 
true  distribution  of  load  in  the  asperity  contact  area.  But  the  one  of 
paraboloidal  distribution  does  give  a  conservectlve  estimate  of  crack 
initiation  and  may  very  well  be  closer  to  the  real  distribution  of  load. 
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Figure  6.  Stress  fields  for  varying  contact  area  and  load  distribution. 


2.2  Contact  Area  Parameters 

The  study  considers  both  the  shape  of  the  contact  area  and'  the  aspect 
ratio  of  the  contact  area.  Since  the  actual  integral  must  be  carried  out 
numerically,  it  is  difficult  to  generate  a  theory  for  all  contact  area 
shapes.  For  practical  purpose  and  computation  feasibility,  the  thermo¬ 
mechanical  principal  stresses  of  a  circular  and  a  square  contact  areas  will 
be  compared.  There  is  very  little  noticeable  difference;  that  is,  the 
differences  are  less  than  the  permissible  numerical  errors.  It  is  therefore 
postulated  that  the  shape  of  the  contact  area  is  not  essential  in  wear 
analysis  of  asperity  excitation. 

The  effect  of  the  aspect  ratio  (t)  of  the  contact  area  covers  the  study 
or  a  square  (or  circular,  t  =  1)  area  ranging  toward  a  line  contact  (t  =  ®), 
where  the  problem  would  be  reduced  to  two-dimensional.  Figure  6  shows  the 
consistent  trend  of  lower  stresses  with  increase  in  aspect  ratios.  To  the 
limit  is  the  two-dimension!  (t  =  «}  solution,  of  which  the  value  is  only 
one-sixth  of  the  three-dimensional  solution  of  stresses. 

It  is  to  be  emphasized  that  the  two-dimensional  and  the  three- 
dimensional  solutions  of  stresses  from  the  mechanical  loading  portion  are 
essentially  the  same.  Hence,  the  phenomenon  is  principally  thermal. 

Indeed,  both  the  two-dimensional  and  the  three-dimensional  theories  predict 
approximately  the  same  surface  temperature.  But  the  three-dimensional  model, 
because  of  the  lateral  (perpendicular  to  the  plane)  heat  loss,  will  drop 
temperature  more  rapidly  thus  higher  temperature  gradient,  leading  toward 
higher  thermal  stresses. 

The  effects  from  different  half-widths  are  shown  in  Figure  7.  It  is 
noted  that  higher  solutions  result  from  larger  asperity  size  which  leads  to 
more  heat  input,  thus  higher  temperature  gradient.  However,  larger  width 
Implies  larger  contact  area  with  correspondingly  reduced  pressure.  The 


Increasing  stress  with  the  Increase  of  half-width  (a)  is  less  than  linear, 
which  may  be  observed  in  Figure  7*,  while  the  decrease  in  pressure  is  most 
likely  inversely  proportional  to  the  square  of  the  characteristic  dimension 
(a).  The  result  is  definitely  beneficial  for  larger  a,  if  it  is  design 
controllable. 


3.0  MATERIAL  PROPERTY  PARAMETERS 


The  material  parameters  under  consideration  are  the  mechanical  consti¬ 
tutive  coefficients  (\,n)  or  (E,v),  respectively  the  Lame's  coefficients  and 
the  pair  of  Young's  Modulus  and  Poisson's  Ratio,  and  the  mass  density  (p), 
the  coefficient  of  Coulomb  friction  (^),  the  coefficient  of  thermal  expan¬ 
sion  (a),  the  thermal  diffusivity  {*•)  and  the  thermal  conductivity  (k). 

Their  effects  are  mostly  thermal.  The  mechanical  portion,  as  a  result  of 
tractions  only,  is  essentially  related  to  the  prescribed  loads.  For  the 
problem  of  the  current  asperity  speed,  it  is  the  thermal  component  that 
dominates  the  thermomechanicaT  cracking  phenomenon.  In  order  to  keep  the 
mathematics  manageable,  some  simplified  postulations  are  necessary.  For 
instance,  the  report  adopts  the  simple  Coulomb's  law  on  the  friction  force, 
Ff  =  nfFn.  The  heating  effect  of  the  friction  force  influences  in  the  sur¬ 
face  layer  a  temperature  rise,  which  In  turn  increases  the  Coulomb  coeffi¬ 
cient  ji-f  with  resulting  further  increase  in  friction  force  [2],  The 
analysis,  therefore,  avoids  the  iteration  by  presuming  a  relatively  high 
constant  coefficient  to  anticipate  a  steady  state  high  temperature  field. 

The  numerical  computation  of  thermal  stresses  is  thus  based  on  a  coefficient 
of  Coulomb  friction  of  0.5,  which  is  conservative  for  high  coefficient  of 
friction.  The  mass  density  (p)  affects  the  mechanical  portion  of  the  stress 
field,  Equation  (18),  through  the  Mach  numbers  (M  )  and  the  thermal  field 
through  the  thermal  diffusivity  (r  =  k/pc).  Since  the  Mach  number  is  gene¬ 
rally  small  in  this  class  of  problems  and  the  mechanical  portion  of  the 
stress  field  is  less  significant,  the  effect  of  the  mass  density  is  princi¬ 
pally  combined  with  the  specific  heat  (c)  to  define  the  thermal  capacity 
(pc)  in  affecting  the  value  of  the  thermal  diffusivity  {<). 


3.1  Mechanical  Properties 

It  is  noticed  in  the  thermal  stress  equation  (21)  and  the  definitions 
of  the  coefficients  (b.)  that  the  thermal  stresses  are  directly  proportional 
to  its  stiffness  (y  or  £)  —  modulus  of  rigidity  or  Young's  modulus,  and  the 
coefficient  of  thermal  expansion  a  (as  a  factor  in  the  coefficient  a^).  The 
effect  of  the  Poisson's  ratio  is  rather  obscure.  Figure  8  illustrates  the 
fact  that  materials  with  five  times  difference  in  stiffness  does  indeed  show 
corresponding  difference  in  the  maximum  principal  stress.  The  variance  of 
the  Poisson's  ratio,  which  is  usually  around  0.3,  does  not  indicate  to  be  a 
significant  factor.  Yet,  as  shown  in  Figure  8,  the  principal  stress  will 
increase  slightly  with  increase  in  the  Poisson's  ratio.  It  is  also  noticed 
that,  even  though  the  mechanical  portion  of  the  stress  components  are 
independent  of  the  material  stiffness.  Equation  (18),  the  effect  of  stiff¬ 
ness  on  the  combined  thermomechanical  principal  stress  is  still  significant, 
owing  to  the  dominance  of  the  thermal  effect  in  high  speed  asperity  exci¬ 
tation. 

3.2  Thermal  Properties 

Similar  to  the  effect  of  stiffness  of  the  wear  material,  the  mechanical 
portion  of  stresses  is  also  independent  of  the  coefficient  of  thermal  expan¬ 
sion.  But  the  combined  thermomechanical  stress  is  almost  directly  propor¬ 
tional  to  the  coefficient  of  thermal  expansion  as  indicated  in  Figure  9. 
Physically  it  is  apparent  that  higher  coefficient  of  thermal  expansion 
introduces  more  mechanical  constraint  in  the  material  for  the  same  temper¬ 
ature  gradient. 

The  effect  of  the  thermal  diffusivity  {<)  is  more  complex.  Although 
the  thermal  diffusivity  occurs  as  the  only  thermal  property  apparent  in  the 
governing  differential  equation  (5),  its  value  is  determined  by  the  coeffi- 
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cient  of  thermal  conductivity  (k)  and  the  thermal  capacity  (pc)  in  that 


r  =  k/pc.  If  the  thermal  conductivity  (k)  is  kept  constant,  the  lowering  of 
thermal  diffusivity  will  lower  the  temperature  field  because  of  increase  in 
the  thermal  capacity  (pc).  It  lowers  the  temperature  field;  and  its  effect 
on  the  thermal  stress,  is  positive.  However,  if  the  reduction  in  the  ther¬ 
mal  diffusivity  is  caused  by  a  reduction  in  thermal  conductivity  (k),  the 
boundary  condition  (Equation  17)  is  also  affected,  resulting  in  higher 
temperature  gradient  thus  higher  thermal  stress.  Figure  10  illustrates  the 
case  for  which  the  thermal  capacity  (pc)  remains  constant.  The  change  in  the 
thermal  diffusivity  is  directly  proportional  to  the  change  in  the  thermal 
conductivity.  The  maximum  principal  stress  increases  with  decrease  of  the 
thermal  conductivity,  thus  thermal  diffusivity,  even  though  not  exactly  in 
inverse  proportion.  However,  if  the  reduction  in  the  thermal  diffusivity  Is 
due  to  the  increase  in  thermal  capacity  rather  than  the  decrease  in  thermal 
conductivity,  lower  stress  will  result.  The  effects  of  variation  in  thermal 
properties  on  the  maximum  thermal  stress  at  the  depth  c  =  0.1  for  a  uni¬ 
formly  loaded  rectangular  contact  area  t  =  10  and  an  asperity  speed  of 
15  m/s,  is  illustrated  in  Figure  11.  Curves  1  through  3  maintain  a  constant 
thermal  capacity.  Curve  1  uses  the  thermal  properties  of  Stellite  III. 

Curves  2  and  3  illustrate  the  decrease  and  increase  of  the  principal  thermal 
stresses  corresponding  to  the  doubling  or  halving  the  thermal  diffusivity. 
However,  when  the  halving  of  the  thermal  diffusivity  is  caused  by  a  doubling 
in  the  thermal  capacity  rather  than  halving  of  the  thermal  conductivity,  the 
principal  thermal  stress  reduces  almost  in  half  as  shown  by  Curve  4. 

Finally,  as  significantly  illustrated  in  Curve  5,  the  thermal  diffusivity  is 
kept  the  same  as  in  Curve  1,  but  both  the  thermal  conductivity  and  the  ther¬ 
mal  capacity  are  halved.  The  ensuing  increase  in  thermal  stress  is  most 
noticeable. 


Effect  of  thermal  conductivity  on  the  stress  field 


1.  x,  =  k  /  CX  =  C 

2.  t<2  *'2k  ,  k2  -  2k  ,  C2  =  c 

3*  x^-x/2  ,  *  k  /  2 03*0 

4.  =  x  /  2  ;  =  k  y  Q|  =  2c 

5,  x^®x<,k5ssk/2>C5  =  c/  2 


Fig.  11.  Effect  of  variation  in  thermal  properties 


4.0  CONCLUSION 


The  parametric  study  enhances  the  understanding  or  thermomechanicai 
cracking  through  the  mechanism  of  subsurface  cohesive  failure,  which  is 
predominantly  a  thermal  stress  phenomenon  at  high  speed  asperity  excitation. 
The  asperity  parameters  are,  in  general,  not  controllable  quantities.  The 
purpose  is  therefore  to  improve  analysis  for  conservative  designs.  For 
instance,  a  parabolic  distribution,  resulting  in  higher  themomechanical 
stress,  is  preferred  over  the  assumption  of  uniform  distribution.  The 
asperity  initiated  cracking  is  essentially  a  three-dimensional  phenomenon; 
two-dimensional  analysis  under-estimated  stress  level  by  six  times.  The  use 
of  uniform  load  distribution  and  two-dimensional  formulation  cannot  be 
justified  for  actual  design  computation  unless  an  artificial  and  high  factor 
of  safety  is  used.  Their  use  may  be  used  in  cases: 

1.  Two-dimensional  theory  yields  the  same  mechanical  stress  field  and 
the  same  surface  temperature  as  those  from  the  three-dimensional  theory 
developed  in  [6,7].  In  such  case,  two-dimension  theory  may  be  preferred  for 
its  much  simpler  mathematical  computation. 

2.  Uniform  distribution  of  load  may  be  employed  for  the  illustration  of 
parametric  effect  in  the  comparison  of  stress  magnitudes.  Since  the  rela¬ 
tive  stress  levels  rather  than  the  real  ones  are  of  interest,  the  simpli¬ 
fication  in  mathematical  analysis  Is  justified. 

The  other  non-control  1  able  asperity  parameter  is  the  asperity  shape.  It 
is  gratifying  with  the  demonstration  of  the  comparison  of  a  square  and  a 
circular  contact  area,  both  of  an  aspect  ratio  of  t  =  1,  that  the  difference 
of  their  results  Is  negligibly  small.  The  asperity  speed  affects  the  thermo- 
mechanical  field  two-fold.  In  one,  the  higher  speed  results  in  a  higher 
heat  input,  VR-, ,  leading  toward  a  higher  thermal  stress  field.  -In  another, 
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the  higher  speed  increases  the  Peclet  number,  causing  a  decrease  of  the 
thermal  layer,  the  distance  between  the  wear  surface  and  the  depth  of  maxi¬ 
mum  temperature  gradient.  This  decrease  in  the  thermal  layer,  itself  would 
cause  an  increase  in  the  temperature  gradient,  thus  promoting  an  increase  in 
the  thermal  stress. 

Among  the  material  properties,  the  stiffness  (n  or  E)  and  the  coeffi¬ 
cient  of  thermal  expansion  (a)  are  directly  proportional  to  the  principal 
thermal  stress.  They  are  unrelated  to  the  mechanical  portion  of  the  stress 
field.  In  view  of  the  predominance  of  thermal  stress  field  for  high  speed 
asperity  excitation,  the  choice  of  stiffer  wear  material  must  be  compensated 
by  a  sufficient  reduction  In  the  coefficient  of  Coulomb  friction. 

The  thermal  properties  to  be  considered  in  design  should  be  the  thermal 
conductivity  (k)  and  the  thermal  capacity  (pc).  The  thermal  diffusivity  is 
a  derivative  parameter.  High  thermal  conductivity  and  high  thermal  capacity 
are  definitely  preferred  to  reduce  the  thermal  stress.  Furthermore,  an 
increase  in  thermal  conductivity  would  lead  to  lower  Peclet  number 
(Pe  =  Vjl/k- )  thus  thicker  thermal  layer.  The  ensuing  reduction  in  temper¬ 
ature  gradient  will  result  in  a  further  decrease  in  the  thermal  stress 
field. 

The  present  report  thus  identifies  those  parameters  to  be  considered  in 
analysis  and  in  design.  The  information  should  assist  the  analyst  to  choose 
the  proper  mathematical  model  and  the  designer  to  select  wear  materials 
based  on  the  Influential  material  properties  for  an  optimal  design. 
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